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I. INTRODUCTION 

Lorentz invariance is one of the fundamental principles of Einstein’s general relativity (GR) and modern 
physics. The success of GR to describe all observed gravitational phenomena, together with its intrinsic 
mathematical elegance is interpreted as a further proof of the importance of the Lorentz invariance [T]. 
However, Lorentz invariance may not be an exact symmetry at all energies [2]. Any effective description 
must break down at a certain cutoff scale signaling the emergence of new physical degrees of freedom beyond 
that scale. For example, the hydrodynamics, Fermi’s theory of beta decay [3] and quantization of GR |4] 
at energies beyond the Planck energy. Astrophysical observations suggest that the high energy cosmic rays 
above the Greisen-Zatsepin-Kuzmin cutoff is a result of the Lorentz violation [5]. Lorentz invariance also leads 
to divergences in quantum field theory which can be cured with a short distance of cutoff that breaks it [5]. 
Often, the late cosmic acceleration is also interpreted as a demand for a modification of GR at cosmological 
scales [71 [5]. 

There are several gravitational theories that violate the Lorentz invariance [2], e.g., Hofava-Lifshitz theory 
0, ghost condensations [TO], warped brane worlds, Einstein-aether theory etc. In Einstein-aether theory, 
the background tensor helds break the Lorentz symmetry, and were once thought must to be dynamical m, 
but more careful investigations recently revealed that it is not necessary [T^. In this theory, the Lorentz 
symmetry is broken only down to a rotation subgroup by the existence of a preferred time direction at every 
point of spacetime, i.e., existing a preferred frame of reference established by aether vector rt“. This timelike 
unit vector held can be interpreted as a velocity four-vector of some medium substratum (aether, vacuum 
or dark huid), bringing into consideration of non-uniformly moving continuous media and their interaction 
with other helds. Meanwhile, this theory can be also considered as a realization of dynamic self-interaction of 
complex systems moving with a spacetime dependant macroscopic velocity. As to an accelerated expansion of 
the universe, this dynamic self-interaction can produce the same cosmological effects as the dark energy m- 

The introduction of the aether vector allows for some novel effects, e.g., matter helds can travel faster than 
the speed of light new gravitational wave polarizations can spread at different speeds |15j . It should 
be noted that the propagation faster than that of the light does not violate causality m- In particular, 
gravitational theories with breaking Lorentz invariance still allow the existence of black holes li inHis]- 
However, instead of Killing horizons, now the boundaries of black holes are hypersurfaces, termed universal 
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horizons mui], which can trap excitations traveling at arbitrarily high velocities. This universal horizon may 
radiate thermally at a fixed temperature and strengthen a possible thermodynamic interpretation though 
there is no universal light cone [20] (See also |2T| for a different suggestion.). 

It is naturally to extend the Einstein-aether theory to include other fields, i.e., electromagnetic one | 22 |- As 
for cosmology, the interaction of electromagnetic waves with a non-uniformly moving aether can change the 
details of the standard history of the relic photons that could be tested using observational data. As for black 
holes, the interaction of electromagnetic radiation with a deformed aether will induce new dynamo-optical 
effects that could be also tested. As for gravitational waves, the Einstein-Maxwell-aether theory is expected 
to predict new forms for gravitational wave propagations mm- Our goal here is to extend Einstein-aether 
theory to include a source-free Maxwell field. For more general formalism of the Einstein-Maxwell-aether 
theory, see [H]. 

The rest of the paper is organized as follows. In Sec. II we provide the background for the Einstein-Maxwell- 
aether theory to be studied in this paper. In Sec. Ill we construct a Smarr formula for spherically symmetric 
solutions. In Sec. IV, we first construct two new classes of exact charged solutions, and then use them as 
examples to study the Smarr formula. In Sec. V, we present our main conclusions. 

Before proceeding further, we would like to note that the exact charged solutions presented in this paper can 
be considered as a generalization of the neutral {Q = 0) ones given in [T^. Therefore, in the following there 
may exist repeating of the materials presented there, in order for the current paper to be as much independent 
as possible, although we shall try to limit this to its minimum. For more detail, we refer readers to |lfi] . 

II. EINSTEIN-MAXWELL-AETHER THEORY 


The general action for the Einstein-aether theory can be constructed by assuming that: (I) it is general 
covariant; and (2) it is a functional of only the spacetime metric gab and a unit timelike vector u°‘, and 
involves no more than two derivatives of them, so that the resulting field equations are second-order differential 
equations of gab and u“. Then, the Einstein-Maxwell-aether theory to be studied in this paper is described 
by the action. 


5 = 




1 


-IGttG: 


-{TZ -t- £ae) + 


M 


In terms of the tensor defined as |23l [24] , 


( 2 . 1 ) 


= Ci5“^5cd + C25\5^^ + - C4u‘^u’’gcd , 


( 2 . 2 ) 
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the aether Lagrangian is given by 

-L^ = - \{v? + 1), (2.3) 

where Ci(i = 1,2,3,4) are coupling constants of the theory. The aether Lagrangian is therefore the sum 
of all possible terms for the aether field u°‘ up to mass dimension two, and the constraint term A(m^ + 1) 
with the Lagrange multiplier A implementing the normalization condition v? = —1. The source-free Maxwell 
Lagrangian Cm is given by 

>Cm = — „ „ ^ab = VaAlb — Vb^Q, (2.4) 

lOTrGa, 

where Aa is the electromagnetic potential four-vector. 

There are a number of theoretical and observational bounds on the coupling constants Ci [3III1I1S]. Here, 
we impose the following constraints^ 

0 ^ Ci4 <2, 2 -|- Ci3 -f 3 c 2 >0, 0 < Ci3 < 1, (2-5) 


where C 14 = ci -I- C 4 , and so on. The constant is related to Newton’s gravitational constant Gjv by 
Gae = (1 — cii/2)GM, which can be obtained by using the weak field/slow-motion limit of the Einstein-aether 
theory [511^. 


The equations of motion, obtained by varying the action (2.11 with respect to gab^ u“, and A are 


Qab = rs, + ^^G^r^, iE„ = o, v“j-,b = o, u^ = -i, 


( 2 . 6 ) 


respectively, where the aether and Maxwell energy-momentum stress tensors and 7/^ are given by 
= AUaMb + C4aaab - -hCi[(VaUc)(VbM'’) - (V“Ma)(VcMb)], 


q-M _ ^ 

47rG„ 


— —n T I TT "C r 

am*' h 


(2.7) 


with 


^a = Vbr\ + An, + C4(V,n'')ab, = + u^Y^ab)- (2.8) 

The acceleration vector a“ appearing in the expression for the aether energy-momentum stress tensor is defined 
as the parallel transport of the aether field along itself, a“ = Vun“, where Vx = Tf^Vb. 


^ Note the slight difference between the constraints imposed here and the ones imposed in m, as in this paper we also require 
that vacuum Cerenkov radiation of gravitons is forbidden m- 
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g°-b = + s“s'’ + g“^ e(b) • e(c) = tibc, 


Following [16) . we first define a set of basis vectors at every point in the spacetime, so that we can project 
out various components of the equations of motion. Let us first take the aether field to be the basis vector. 
Then, pick up two spacelike unit vectors, denoted, respectively, by and n“, both of which are normalized 
to unity, mutually orthogonal, and lie on the tangent plane of the two-spheres B that foliate the hypersurface 
Tijj. Finally, let us pick up s“, a spacelike unit vector that is orthogonal to u“, m“, n“, and points “outwards” 
along a S[/ hypersurface, so we have the four tetrad, = (u“, s“, m“, n“), with 

(2.9) 

where -1- n“n*’. By spherical symmetry, any physical vector has at most two non-vanishing 

components along, respectively, and s“, i.e., = Aiu‘^ + A 2 s‘^. In particular, the acceleration a“ has only 

one component along s“, namely, a°' = {a ■ s)s“. Similarly, any rank-two tensor Fab may have components 
along the directions of the bi-vectors UaUb, U(^aSb), U[aSb], SaSb, gab, where gab is the projection tensor onto 
the two-sphere B, bounding a section of a hypersurface. In the following, we study the expansion of 
the Maxwell field Killing vector surface gravity k, energy-momentum stress tensors and , 
and Ricci tensor TZab- The given source-free Maxwell field can be formulated in terms of four-vectors 
representing physical fields. They are the electric field E°‘ and magnetic excitation 5“ as. 


^abmn 


= F'^^ub, = 




Fmn'^b, 


( 2 . 10 ) 


where is the Levi-Civita tensor. From Eq.(2.6l it can be shown = 0. Then, we find 


F^b ^-E'^u^ + E^u‘^. (2.11) 

On the other hand, the electric field is spacelike, since E'^Ua = 0. So, we have = {E ■ s)s“. Thus, 


Fab = {E ■ s){—SaUb + SbUa)- After substituting it into (2.6), we can see {E ■ s) = where Q is an integral 


constant, representing the total charge of the space-time. Therefore, we have 

T -Qf _ t 
ab — ^2 \^o.^b '^b^aj- 


( 2 . 12 ) 


The Einstein, aether and Maxwell equations of motion (2.6) can be decomposed by using the tetrad defined 
above. In particular, the aether and electromagnetic energy-momentum stress tensors and the Ricci tensor 
can be cast, respectively, in the forms, 

Eab ~ 'T'au'^aUb — ‘^FasU{aSb) + EggSaSb H ^gab, 

F-ab F-au^a^b FEas'^i^a^b') “b F-ggSaSb T ^gab, 


Tm 


7 -M n-M I n-'M , > ivi ^ 

ab — 'uu'^aUb — ^lus '^{aSb) + Os ^aSb H 'i^gab- 


(2.13) 
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The coefficients of Tjf and in (2.131 can be computed from the general expression (2.7). The corresponding 
coefficients for TZab, on the other hand, are computed from the definition [Va, t)y choosing 

or s“, and then contracting the resulting expressions again with u°‘ and/or s“ appropriately. The 
coefficients for the three (u, s) cross terms are 


r“=ci 4 i^(a-s)+v„(a-s) , T” = 0, TZ^s = (Kq - K/2)k - VsK, 


(2.14) 


where 


^[a^6] — k — g /ls§ab: ^ — ^9 ^uOabi 


(2.15) 


2 ^ —*^/u,L /5 -- — 2 *- 

with K (= Kq + K) being the trace of the extrinsic curvature of the hypersurface S[/. The aether equation 
s • jE = 0 and the us-component TZus = T^s + 87rGae7)(^ yield 


cmVsXo - (1 - ci3)(ifo - K/2)k + (1 + C2)V,i7 = 0, 

Cl23V,X-(l-Ci3)7;*=0. 


(2.16) 

(2.17) 


The uu- and ss-components of the gravitational field equations give 


/-t)l 


(fc + Vs)(a • s) + — (1 — C 13 ) I Kq H —— 


C 123 J.2 


- (1 + C2 + V„K - ^K^ - ^ = 0, 


2 / 


C123 


{K + Vu){,K — Kq) + (1 + C 2 )KKq ~ (1 4—[Vs(a • s) + a^] 


0 k Ci4 / \ 

-\-c14a — Vs 2 ^ ^ 




In the next sections, we will use these equations to obtain new black holes solutions. 


(2.18) 


(2.19) 


III. SMARR FORMULA 


The studies of black holes have been one of the main objects both theoretically and observationally over 
the last half of century [351 US] j and so far there are many solid observational evidences for their existence in 
our universe. Theoretically, such investigations have been playing a fundamental role in the understanding 
of the nature of gravity in general, and quantum gravity in particular. They started with the discovery of 
the laws of black hole mechanics m and Hawking radiation |31j . and led to the profound recognition of 
the thermodynamic interpretation of the four laws |32] and the reconstruction of general relativity (GR) as 
the thermodynamic limit of a more fundamental theory of gravity [33j . More recently, they are essential in 
understanding the AdS/CFT correspondence [34l|35] and firewalls [36| . 
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To derive the Smarr formula, we first introduce the ADM mass, which is identical to the Komar mass 
defined in stationary spacetimes with the time translation Killing vector m, 

Madm = -y^ f ^^X^’d^ab, (3.1) 

dTrGae Jis^ 

where d^ab = ~U[aSb]dA, with dA (= smOdOdip) being the differential area element on the two-sphere B, 
and Boo is the sphere at infinity. The derivative of the Killing vector x“ = — (u • + (s ’ x)®^ is given by 

V“x^ =(3.2) 


where n denotes the surface gravity usually defined in GR, and is given by 


^ = \/“2 


(3.3) 


At the infinity, we have {u ■ x) = —1 and (s • x) = 0- Then, Eq.(3.1) yields. 


Madm = lim 

r—^oo 


r'^{a ■ s) 
Gee 


(3.4) 


In the studies of black hole physics, the physics of horizons has provided useful information. In particular, 
at the Killing horizon the first law m and Smarr formula [3H1 for the Reissner-Nordstrom black hole take the 
forms, 

hkhSAkh , VkhSQ kkhAkh , VkhQ 

oMadm = —-H--, Madm = . ^ -f -, (3.5) 

oTTLzjv Crjv 47r(_TjV 

where Madm is the ADM mass of the spacetime and, kkh[= «^(?’iCir)], Akh and Vkh are the surface gravity, 
cross-sectional area and electromagnetic potential evaluated on the Killing horizon, respectively. Identifying 
Tkh = kkh/‘2'!t as the temperature of the horizon and, the entropy S = Akh/^Gn-, one can obtain the 
analogy with the first law of thermodynamics, SE = T6S + V6Q and E = 2TS + VQ. 

Any causal boundary in a gravitational theory should have an entropy associated with it. Therefore, in the 
Einstein-aether theory, the universal horizons are expected to have also their entropy and the first law of black 
hole mechanics, though whether this entropy is to be proportional to its area or not is still an open issue. 
Meanwhile, these black holes still have Killing horizons. Then, a question is: How is their thermodynamics? 

To obtain some hints, in this section we shall present the Smarr formulas of the universal and Killing horizons 
for general static and spherically symmetric Einstein-Maxwell-aether black holes. Let us first consider the 
geometric identity [30] . 


7^a6x'' = V'’(VaXfc). 


(3.6) 










From the Einstein field equations (2.6), we find that TZac = — gabT’^/2 + 8 ttGs.TJ^, where [3] 


SnG^Tj^X^ = -%xa, 

XS, - x^ = V'’|[ci4 (a ■ s){u-x) + {C 123 K - 2 C 13 K 0 ) (a ■ s)]'U[aS6]|. 


(3.7) 


Setting 


- ^Xa = 2V'’ (E<3(r)n[„S6]) , 


(3.8) 


we find that Eq.(3.6) can be cast in the form, 

= 0, F“'' = 2F(r)u[“s'’l, 

where 


(3.9) 


F(r) = F^{r) + q{r), 
q{r) = - (1 - 


C 123 


2 y (« • Xiu ■ X) + [(1 - £13)7^0 + -^K (s • x). 


(3.10) 


On the other hand, comparing Eq.(3.9) with the soucre-free Maxwell equations (2.6), we find that its solution 
must also take the form (2.12), that is, F{r) = F^jr^. To determine the integration constant Fq, we note 


that, for asymptotically flat space-time, we have [ima, 


ro 


M-X^-1, s-x-0, a-s=^+ 0 {r^), 


(3.11) 


as r —>• 00 . Then, from Eq.(3.10) we find that Fq = ro(l — ci4/2)/2. Thus, we have 


(3.12) 


Inserting Eq.(3.11) into Eq.(3.4), on the other hand, we find that the ADM mass is given by Madm = 


ro/2Gae- Therefore, the total mass M of the spacetime is 

M . + M. = (1 - f) ^ = F4A. 


(3.13) 


where Mg. = —ci/^Madm/‘^ is the aether mass or aether contribution to the renormalization of Madm 


On the other hand, using Gauss’ law, from Eq.(3.9) we find that 


0= f {VbF‘^'^)dEa= f F‘^^di:ab- f F<^'^d^ab= f FdA- f FdA. (3.14) 

J'S JBao JBh dBac JBh 

Here is the surface element of a spacelike hypersurface E. The boundary 9E of E consists of the boundary 


at spatial infinity Boo, and the horizon Bh, either the Killing or the universal. Note that Eq.(3.14) is nothing 
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but the conservation law of the flux of Substituting the above expression into Eq.(3.13) and taking 


Eq.(3.10) into account, we And the following Smarr formula, 


MG. = + v„„Q. MG. = + VkhQ. 


( 3 , 15 ) 


Att Att 

where Ajjh and Akh are, respectively, the area of the universal and Killing horizons, and M is the total mass 
of an asymptotically flat solution defined in the asymptotic aether rest frame. The potential Vh is defined 
as Vh = Hence, the first law for the aether black hole may be obtained via a variation of 

the Smarr relation. In the next section we consider it for two new classes of exact charged aether black hole 
solutions. 

Eor the surface gravity at the universal horizon, when one considers the peeling behavior of particles moving 
at any speed, i.e., capturing the role of the aether in the propagation of the physical rays, one finds that the 
surface gravity at the universal horizon is [MUSI HZ] 


KUH = ■x)= ^ (a • s) (s • x) 


(3.16) 


where in the last step we used the fact that Xa is a Killing vector, V(aXb) = 0- H niust be noted that this is 


different from the surface gravity defined in GR by Eq.(3.3). In particular, at the universal horizon we have 


M • X = 0, and Eq.(3.3) yields 


K(ruH) = F^ois- X)\r=' 


r^rjjH ’ 


(3.17) 


IV. EXACT SOLUTIONS OF CHARGED AETHER BLACK HOLES 

To construct exact solutions of charged aether black holes, let us first choose the Eddington-Einklestein 
coordinate system, in which the metric takes the form 


ds^ = —e{r)dv^ + 2f{r)dvdr + 


and the corresponding timelike Killing and aether vectors are 


(4.1) 


X“ = (l,0,0,0), = (a, /?, 0, 0,), = (-ea + //3, /a, 0, 0,) I I , (4.2) 

d4>. 

where a{r) and /3(r) are functions of r only. Then, the metric can be written as gab = —UaUf, + SaSb + gab, 
where we have the constraints = —1, s^ = l, u - s = 0. The boundary conditions on the metric coefficients 
are such that the solution is asymptotically flat, while those for the aether components are such that 


lim = {1, 0, 0,0} . 

1—yoo 


(4.3) 
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Some quantities that explicitly appear in Eqs.(|2.16l-(2.19| are 


^ ^ {u- xY {s ■ xY 

[a-s) = -^—, Ko = -- 


K = - 


2 (s- x) 


k = - 


2 (m • x) 


/ ’ / ’ “ r/ ’ •“ rf 

where a prime (') denotes a derivative with respect to r. And a{r)^ j5{r) and e(r) are 

1 


a{r) = 


(s • X) - (w • x) 


, P{r) = e{T) = (u-xf-{s-xf- 


Then, from Eqs.(3.16) and (3.17) we obtain 


1 


i^UH = - y{u-xY {s- x) 


f 


, n{ruH) = - 7 (s • X)' (s • X) 

UH J 


(4.4) 


(4.5) 


(4.6) 


UH 


Clearly, in general kuh u-ifun)- 

Erom the above expressions one can see that all quantities can be calculated from (tt • x) and (s • x) under 


(4.7) 


the condition /(r) = 1. A straightforward calculation of Eq.(2.14) yields 

_ 2(s-x)(m-x)/'(0 

rfHr) 

In the static spherical symmetric and asymptotically flat spacetime, if we assume that / = 1 holds in the 
whole space-time, we find 


T^us = TZ = 0, if = !)• 


(4.8) 


From Eq.(3.8), we also find that 


= = 5 . (/ = !)■ («) 

In the following, we shall use the above expressions first to obtain two classes of exact solutions for the 
cases Ci 4 = 0, C 123 Y" 0 and C 123 = 0, C 14 Y 0; all with / = 1. Then, we shall study their main properties by 
using the Smarr formulas given above. 


A. Exact solutions for C14 = 0 


When the coupling constant C 14 is set to zero and C 123 Y Oj from Eqs.(2.17) and (4.8) one can see the 


quantity VsAl has to be vanished, i.e., VgAT = 0. So, the trace of the extrinsic curvature K of the Sy 
hypersurface is constant. In the infinity, this constant will vanish asymptotically due to the asymptotical flat 


conditions. Therefore, it must vanish everywhere. Substituting K = 0 into Eqs.(2.18) and (2.19), we obtain 


Y-x) = ^, 




( 4 , 10 ) 
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where rg, is another integral constant. Then, using the formula (4.5), we find 


7* 7*^ 7*^ 


(4.11) 


which reduces to those given in Ref. [16] when Q = 0. 

The location of the universal horizon rjjH is the largest root of equation u • % = 0. Meanwhile, u • x is a 


physical component of the aether, and should be regular and real everywhere. However, from Eq.(4.10l one 


can see that in the region r_ < r < run, this term becomes purely imaginary, where r_ is another root of 
?x • X = 0, unless the two real roots coincide. Then, rg, becomes a function of tq. That is, the global existence 
of the aether reduces the number of three independent constants (rg, rg,, Q) to two, (ro,(3), the same as in 
GR. Thus, from {u ■ x)^ = 0 and d{u ■ xY I dr = 0 [TH], we find 


ro /3 

= T 4 




'''UH ~ ’ 


(4.12) 


, ,16 1 - ci3 

which is showed in Figure One can see that the charge Q is subjected to the condition Q < 3ro/4-\/2, in 


''uh/^'o 


^kh/^0 




Q/ro 


FIG. 1. The universal and Killing horizons of the charged aether black hole with different C 13 in the case of C 14 = 
0, C123 7^ 0. The presence of the charge Q makes both horizons smaller. The universal horizon does not depend on 
C13, while the Killing horizon becomes bigger with the increasing of C13. When C13 = 0, the Killing horizon reduces to 
that of the Reissner-Nordstrom black hole. 


order to have tuh real. When Q = ro/2, we find rg, = 0 and rjjH = tkh = ?'o/2. When Q > ro/2, we have 
< 0. Thus, in order to have the aether be regular everywhere, the charge should be. 


Q < 


which is the same as that given in the Reissner-Nordstrom black hole. 


(4.13) 


Now let us derive the Smarr relation. Using Eq.(3.16), the surface gravity at the universal horizon can be 
computed and is given by 

1 / 




1 - 




UH 


1 - 




' UH 


(4.14) 
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which is showed in Figure When Q = 0, we find that run = 3ro/4 and kjjh = 3 ( 111 ^ 7^1 which is the 

same as those given in uniiizi- The Smarr formula at the universal horizon is 


'^uh/^o 



0.0 0.1 0.2 0.3 0.4 0.5 


^■kh/^o 



Q/ro 


FIG. 2. The surface gravity at the universal and Killing horizons of the charged aether black hole in the case ci 4 = 
0, C123 7^ 0. To compare with the results given in Ref. m, we shift the kuh line with C13 = 0.1 from the left-hand 
side as a dot-dashed line to the right. One can see that when C 13 is small, kuh is larger than kkh in the low charge 
region, while lower in the large charge region, similar to that given in Ref. m for the Khronometric theory. 


MGsi = — — - ^VuhQi < 1 uh = - 73 — ^uh = — 


-ruH r 


UH 


ruH 


(4.15) 


which don’t depend explicitly on the coupling constants c^’s, because now we have C 14 = 0 and = 0 . It 


is easy to see that quh isn’t proportional to kjjh given by Eq.(4.14). On the other hand, at the universal 
horizon we find 


G^5M = T / 2 


Stt V iruH 5rfj 


r)2 \ 2 

' 6Auh + t,^uh5Q. 


UH 


(4.16) 


Why is there the factor 2/3 in the front of Vuh^ For a better understanding, let us use the method proposed 
in m. i e., using M’s expression from A = 


1 A Att 

G.M=-^-+4Q^ + -Q4, 


(4.17) 


and writing the variation of M as Gg^SM = T5A -\- VSQ, we obtain, 


T = 


d{Gg,M) 1 / 2 




dA 


Stt V 3ruH 3r/. 


UH 


^^_d{G^M) 2Q 2 


(4.18) 


which are the same as those given in Eq.(4.16). However, such defined temperature T is also not proportional 


to Kuh given by Eq.(4.14|. 


On the other hand, the location of the Killing horizon is the largest root of e(r) = 0. Using Eq.(4.11), we 
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find 




L = 




3r" 


21/3(12/+ QV^4) 

3i? 


H 


I = -- 


Cl3 


3-21/3’ - 1-C13 

.X 1/3 


UH 


' UH 

2r3 


H=(^J+ ^-4(121+ Qyrl)3 + J^'^ , J = 271- 721Q^{rl + 2Q'^jrl, 


(4.19) 


which is showed in Figure When C 13 = 0, we find that tkh = ^+j that is, it coincides with the Reissner- 
Nordstrom black hole Killing horizon (here after we denote r± = (1 ± \/7 — 40^lr^rQl2.). Then, the Smarr 
formula and surface gravity (using Eq.( |3.3| )) at the Killing horizon are 

2 3ro , 


(IkhAkh , Tr ^ I'D 

-MG* = -^-h VkhQ, Qkh = 


, KkH = 


r.3 ’ 


(4.20) 


47r ' V 2r2 2r2 

KH 'kH/ ' KH KH ' KH 

which is showed in Figure Note again that the qkh is still not proportional to the kkh- When C 13 = 0, 
both qkh and kkh reduce to those given in the Reissner-Nordstrom black hole, qkh = kkh = (r'+ —?"-)/2r^. 
The first law at the Killing horizon cannot be obtained via the variation method, although it may be obtained 
via Smarr’s method m- However, due to its complexity, we shall not consider this possibility, as even we do 
it, we do not expect to get much from such complicated expressions. 


Note that, when C 13 1, from Eq.(4.11) we find that the solution reduces to the usual Reissner-Nordstrom 


black hole with a universal horizon given by (4.12) that is always inside its Killing horizon teh = which 


is the same as that derived in the Khronometric theory [T5| . 

Einally, let us turn to Figure from which we can see that the presence of the charge Q always makes the 
surface gravity kuh lower, while the presence of the constant C 13 always makes it bigger, after the constraints 


(2.5) are taken into account. For the kkh, the situation becomes more complicated. In particular, when both 


C 13 and g are small, the effects of them is similar to that presented in kjjh as shown in the figure. But for 
large C 13 , e.g. C 13 = 0.92, 0.975, the presence of the charge increases the temperature at the beginning and 
then decreases it when the charge becomes very large. 

B. Exact solutions for C 123 = 0 


In this case, from Eq.(3.13l we find that the total mass is 

"°- = (‘-t)t 


(4.21) 


There exists a range of the coupling constants that passes all the current observational tests in the one- 


parameter family of the Einstein-aether theories [40]. Setting C 123 = 0, from Eqs.(2.18), (2.19) and (4.4) we 
























obtain 


14 


(u-x) = -l + g, 


(s • X) = 


To + 2r^ 
2r 


ro 


2 — Ci4 


4Q2 


2(1 - C 13 ) (1 - Ci3)r, 


-1 . 


Then, we find that 


/ X ro r„(ro + r„) 

e(r) = l- —--, /(r) = 1, 

which again reduces to that given in Ref. [16j when Q = 0. From the above expressions, we find 

a{r) = 


(4.22) 


(4.23) 


(4.24) 


(s • X) - (u • x) 1 + V ' 

Since it is one of the component of it should be regular everywhere (possibly except at the singular point 
r = 0), we must have 

^ ^ n ^ /2 ci3 - C14 tq f cu 7 ( ro C14 

r„>0^Q<W- - - — = - — - (1 - C 13 ) y, ci 3 >—. (4.25) 


The position of the universal horizon run and its surface gravity (using Eq. 3.16) are 
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ro 

ruH- rcuH - ' 


7 _ (ill _ 


UH 


And the kuh is showed in Figure]^ Also, when Q = 0 it reduces to the one obtained in uniiiT]. One 


(4.26) 


can see 
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FIG. 3. The surface gravity at the universal horizon and Killing horizon of charged Einstein-aether black holes in the 
case of ci4 = 0.1, C123 = 0. The restrictions are C13 > 0.05, Q/ro < ■y/(2ci3 — 0.1)/8. 


that rjjH does not depend on the charge Q, but kuh depends on it. In other words, to the same universal 
horizon rjjH = ?'o/2, there are different thermal temperatures of the horizon with different charge Q, if we 
assume that T is still somehow proportional to kuh- 
The Smarr formula new reads, 


^ . qUH^UH , T, ^ 1 A Ci4 

GisM = — -h VuhQ, quH = - 1 —s-^ 

47 r ruH \ 2 rf 


UH 


V 2 / ruH oTT 


(4.27) 
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in which the term proportional to SQ is absent. To see this more clearly, let us consider the Smarr method, 
from which we find that 


G^M = 



(4.28) 


Then, writing the variation of M as Gg^SM = T5A + V5Q, we obtain, 

d{Gg,M) 


^_ 9(G3eM) _/ Ci4\ 1 

dA V 2 / Sirru h ’ 


dQ 


= 0 , 


(4.29) 


which are the same as those given in the second line of Eq.(4.27). Once again, qun and T aren’t proportional 


to the KuH given by Eq.(4.26). 


On the other hand, the Killing horizon and its surface gravity (using Eq. 3.31 are 

2ru + ro 


fKH = roAru, kkh = 


KH 


(4.30) 


which are showed in Figure|^ In order for them to be real we must assume that Q < -^1 — Comparing 

Eq.(4.25) with ( 4.30[ ), one can see that the latter condition on Q is contained in the former, i.e., in the charged 
aether black hole, the electric charge is subjected to more stringent restrictions. Again when C 13 = C 14 = 0, 
kkh reduces to (r+ — r_)/2r^. The Smarr formula at the Killing horizon is 


G.M = ,™ = 

Att 


\( 

1 

r-o 

1 

[v 

2 J 

KH 

^3 

' KHJ 


(4.31) 


which depends on the coupling constants C 13 and C 14 . qKH approaches to (r_|_ — r_)/ 2 r^, if C 13 = C 14 = 0 . 
We find that taking variation with respect to each term cannot obtain the first law, so instead we use Smarr’s 

method 133 , and hnd that 


A , DM CaCbQ 

Gg,6M = —6 Akh + — = =, 

'9^ ■\/co(caC6 - l)Q^ + 

dM 


1 - Ci3 


^ _ Cu rp _ _ Cb ^ __CoCb___ ^ 

2’ dA CaCb - 1 \^/ca{CaCb - l)Q'^ + CaCtrj^^ rKHj 


(4.32) 


Note that, similar to the previous case, now qxH and T aren’t proportional to the kkh given by Eq.(4.30), 
either. 


Finally, from Figurej^we note that the dependence of kjjh on Q is similar to the former case. In particular, 
its presence always makes the temperature lower, while the presence of C 13 increases it. The surface gravity 
Huh is always larger than kkh- At the Killing horizon, the effects of the charge and C 13 on khh are just 


opposite. 
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V. CONCLUSIONS 


In this paper, we have studied the Einstein-Maxwell-aether theory, and found two new classes of charged 
black hole solutions for the special choices of the coupling constants: ( 1 ) C 14 = 0 , C 123 ^ 0 , and ( 2 ) C 14 ^ 
0: C 123 = 0. In the hrst case, the universal horizon depends on its electric charge Q, while it doesn’t in the 
second case. In both of the cases, the universal horizons are independent of the coupling constants q, while 
the Killing horizons depend on them. When C 13 (= ci + C 3 ) is very small and approaches to zero, the solutions 
in the case C 14 = 0, C 123 ^ 0 reduce to the usual Reissner-Nordstrom black hole solution. The corresponding 
properties at the universal horizons are the same as those presented in |18j via Khronometric theory. 


To study the solutions further, we have considered their surface gravity and constructed the Smarr formula 
at each of the horizons, universal and Killing. We have shown that there is no problem for such constructions, 
but when trying to construct the corresponding first law of black hole mechanics, they are all different from 
the usual one. In particular, we have shown that the temperature obtained from the Smarr mass-area relation 
is not proportional to its corresponding surface gravity, when both of the charge and aether are present, in 
contrast to the case without aether {a = 0 ) [TH], or to the case without charge m- In particular, in m 
it was found that in the neutral case qjjH is always proportional to the surface gravity kuh at the universal 


horizons, even when the aether is present. In the present paper, from Eqs.(4.201, (4.30) and (4.31) we can 
see that, when Q = 0, is also proportional to hkh- Then, one can also be able to construct a slightly 
modified Hrst law of black hole mechanics at the Killing horizons. 


However, when the charge Q is different from zero, comparing (quH, qKu) with (kuh, kkh), one can see that 
these proportional relations don’t hold any longer. It was also found that in the presence of the cosmological 
constant [31], quH is also not proportional to kjjh- Therefore, it is not clear how to build the first law for 
these charged aether black holes, before we have a better understanding of the entropy of the universal and/or 
Killing horizons. 

The solutions presented in this paper can be generalized to the case coupled with the cosmological constant 
A, which are given by 


e(r) = 


I - 


rols+2r^ 


K 


+ ^-ci3^ + (222^-iA)r2 


I_ ro + 


C 14 — 2 ci 3 '^q 

2(1 —C13) 4 r 2 


1-C13 y 3 / I r ’ 


(ci4 = 0, C 123 yf 0), 
(ci 4 yf 0, C123 = 0), 


(5.1) 


where tq, 1^ and are integration constants. 


In addition, from these solutions, one can also construct topological charged Einstein-aether (anti) de Sitter 
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black holes, which are 


ds^ = — [e(r) — l] dv^ + 2dvdr + (^d9^ + d(j)^'j , 


(5.2) 


where e(r) is given by Eq.(5.1). The studies of the properties of the above solutions are out of the scope of 


this paper, and we hope to report them in another occasion soon. 
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